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ABSTRACT 

The dense, cold gas of Infrared Dark Clouds (IRDCs) is thought to be representative of the initial 
conditions of massive star and star cluster formation. We analyze ^■^CO J = 1 — line emission data 
from the Galactic Ring Survey of Jackson et al. for two filamentary IRDCs, comparing the mass surface 
densities derived from ^'^CO, E13CO1 with those derived from mid-infrared small median filter extinction 
mapping, Ssmf, by Butler & Tan. After accounting for molecular envelopes around the filaments, we 
find approximately linear relations between E13C0 and Ssmf, i-e. an approximately constant ratio 
Si3Co/5^SMF in the clouds. There is a variation of about a factor of two between the two clouds. We 
find evidence for a modest decrease of Siaco /^smf with increasing E, which may be due to a systematic 
decrease in temperature, increase in importance of high -'^'^CO opacity cores, increase in dust opacity, or 
decrease in ^"^CO abundance due to depletion in regions of higher column density. We perform ellipsoidal 
and filamentary virial analyses of the clouds, finding that the surface pressure terms are dynamically 
important and that globally the filaments may not yet have reached virial equilibrium. Some local regions 
along the filaments appear to be close to virial equilibrium, although still with dynamically important 
surface pressures, and these appear to be sites where star formation is most active. 

Subject headings: ISM: clouds, dust, extinction — stars: formation 

1. INTRODUCTION 

Massive, high column density Infrared Dark Clouds (IRDCs), typically identified as being opaque against the Galactic 
background at ^ 10/im, are thought to contain the sites of future massive star and star cluster formation (e.g. Rathborne 
et al. 2006), since their densities (nn > 10* cm""^) and mass surface densities (S > 0.1 gcm~^) are similar to regions 
known to be undergoing such formation activity (Tan 2007). Studies of molecular line emission from IRDCs can help 
determine their kinematics. In particular, we would like to know if they are gravitationally bound, if they are near virial 
equilibrium and if there is evidence for coherent gas motions that might indicate that IRDC formation involves converging 
atomic flows (Heitsch et al. 2008) or converging molecular flows from cloud collisions (Tan 2000). 

In this study we use ^'^CO J = 1 — line emission data from the Galactic Ring Survey (GRS) (Jackson et al. 2006) for 
two filamentary IRDCs, clouds F (^ = 34.437°, h = 0.245°, d = 3.7 kpc) and H (Z = 35.395°, b = -0.336°, d = 2.9 kpc) 
from the sample of 10 relatively nearby massive and dense IRDCs of Butler & Tan (2009, hereafter BT09), comparing 
^•^CO-derived mass surface densities, E13CO1 with small median filter (SMF) mid-infrared (MIR) (8/im) extinction mapping 
■ derived mass surface densities, ESMF7 using the method of BT09 applied to the Spitzer Infrared Array Camera (IRAC) 
' band 4 images of the Galactic plane taken as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) 
. (Benjamin et al. 2003). We consider systematic errors in each of these methods, which is necessary before analyzing larger 
samples of clouds. We are also able to look for evidence of changing CO abundance with column density, e.g. due to 
possible depletion of CO at high densities. We then perform a virial analysis of the clouds to determine their dynamical 
state. 

There have been a number of other studies comparing ^'^CO derived mass surface densities with those from other 
methods. For example, Goodman, Pineda & Schnee (2009) compared near infrared (NIR) dust extinction, far infrared 
(FIR) dust emission and ^^CO line emission in the Perseus giant molecular cloud (CMC), probing values of E up to 
^ 0.02g cm~^ (i.e. up to Ay — 8 mag). Even after accounting for temperature and optical depth variations they concluded 
that ^■^CO emission was a relatively unreliable tracer of mass surface density, perhaps due to threshold, depletion and 
opacity effects. Our study probes higher values of E, from ^ 0.01 to ~ 0.05 gcm~^, and compares ^'^CO emission with 
MIR extinction in order to investigate these processes. 

Battersby et al. (2010) used ^"^CO emission, MIR extinction and FIR dust emission methods to measure E and mass of 
clumps in 8 IRDCs, one of which is IRDC F of our study. They did not present a specific comparison of E13C0 with other 
methods, although derived clump masses were in reasonable agreement. Their sample also included MIR-bright regions, 
associated with ultra-compact HII regions, for which the MIR extinction method cannot be applied. As we describe 
below, our approach differs in a number of ways, including by focusing on filamentary and mostly quiescent regions of 
IRDCs for which the MIR extinction method is most reliable and which are likely to be closer to the initial conditions of 
the massive star and star cluster formation process. We note that while IRDC F in particular does contain some regions 
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Fig. 1. — (a) Left panel; Dependence of with Tb^u (eq. |2} assuming Tex = 10,15,20 K (red dotted, blue solid, green dashed lines, 
respectively), (b) Right panel: Dependence of dNi^cQ /dv with Tg ,^ assuming Tex = 10, 15, 20 K (red dotted, blue solid, green dashed lines, 
respectively) . 

of quite active star formation, including an ultra-compact H II region, here we have concentrated on its more quiescent 
portions. 



We evaluate the column density of ^^CO molecules, dA^iaco: in a velocity interval dv from their J = 1 ^ emission 
via 



where Qrot is the partition function, A = 6.3355 x 10~^s~^ is the Einstein coefficient, Aq = 0.27204 cm, gi = 1 and 
(7„ = 3 are the statistical weights of the lower and upper levels, Ti, is the optical depth of the line at frequency v, i.e. 
at velocity v, T^x is the excitation temperature (assumed to be the same for all rotational levels). For linear molecules, 
the partition function is Q^ot = SjLo(^'^~'~ l)6xp(— i?j/A:Tcx) with Ej = J{J + l)hB where J is the rotational quantum 
number and B = 5.5101 x 10^° s'^ is the rotational constant. Thus for "CO(l-O) we have Ej/k = 5.289 K. For J = 1, 
Qrot = 4.134, 6.018, 7.906 for Tex = 10, 15, 20 K. 
The optical depth is determined via 



where Tb^i, is the brightness temperature at frequency ly, f{T) = [exp(/ii'/[fcT]) — 1]""'^, and Tbg = 2.725 K is the 
background temperature. Tb.v is derived from the antenna temperature, Ta, via Ta = fy/ciump?B,!/, where rj is the main 
beam efficiency (77 = 0.48 for the GRS) and /dump is the beam dilution factor of the ^^CO emitting gas, which we assume 
to be 1 for the IRDCs we are studying. However, it should be noted that the BT09 MIR extinction maps of the IRDCs 
do show that there is structure on scales smaller than the angular resolution of the GRS survey. Given the observed Tb,v 
and for an assumed Tex we solve equations ^ and ([2]) for and thus dA^i3co/df (see Fig. [1]). For a given Tb,,^, this 
figure allows us to judge the sensitivity of the derived column density per unit velocity to temperature uncertainties. 

While we use the above formulae, which account for optical depth, to calculate our column density estimates, for 
convenience we also state their behavior in the limit of optically thin conditions. We have Tb.v — ihiy/k)[f{Tcx)~ f{Thg)]Ti, 
so that 
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The last coefficient changes to (0.9865, 1.347) x lO^^ for Tc^ = 10, 20 K. 

Devine (2009) estimates a temperature of 19 K for cloud F based on VLA observations of NH3(1, 1) and (2, 2). For this 
cloud, we thus adopt a temperature of 20 K. For cloud H we choose a more typical IRDC temperature of 15 K (Carey et 
al. 1998; Carey et al. 2000; Filial et al. 2006). We also consider the effect of varying these temperatures. 

The GRS has an angular resolution of 46", with sampling every 22". The velocity resolution is 0.22 kms~^ (Jackson et 
al. 2006). From a morphological examination of the ^'^CO emission in I, b, v space and comparison to the MIR extinction 
maps of BT09 we identify the velocity range of the gas associated with each IRDC filament (see Figs. [H [3]and|3]). For 
cloud F we consider associated gas to be at LSR velocities 48 — 65 km s~ and for cloud H at 40 — 50 km s^ . The total 
^■^CO column is then evaluated over the velocity range of the cloud iViaco = / dNi^co- 

To convert from A^isco to total mass surface density Sisco we first assume 

1^^6.2:^ + 18.7 (5) 
ni3CO kpc 

(Milam et al. 2005), where i?gai is the galactocentric radius. For clouds F and H we estimate i?gai = 5.37,5.89 kpc, 
respectively (assuming i?gai,o = 8.0 kpc), which would yield ni2co/"'i3CO — 52,55, respectively. For simplicity we adopt 
'^■i2CO /"13CO = 54 for both clouds. We then assume 

= 2.0 X 10-^ (6) 

similar to the results of Lacy et al. (1994). The observed variation in this abundance in GMCs is about a factor of two 
(Pineda, Caselli, & Goodman 2008), and this uncertainty is a major contributor to the overall systematic uncertainty in 
our estimate of Siaco- Thus our assumed abundance of ^^CO to II2 is 3.70 x 10^^ and 

Ei3co = 1.26 x 10-2 J^H2%^g cm-^ (7) 
10^°cm 

assuming a mass per H nucleus of /in = 2.34 x 10^^^ g. 

The assumption of LTE breaks down if the density of the gas is lower than the effective critical density, ^ricrit, i-e. the 
critical density of ^'^CO (J=l-0), ricrit — 1-9 x lO'^ cm~"^, allowing for radiative trapping with escape probability /? — e^^, 
where we approximate f as the column density weighted mean value of t^,. Note, for a spherical cloud we would set f equal 
to the average of but we use for these clouds with a filamentary geometry since we expect t^, seen along our line of 
sight to be relatively small compared to other viewing angles. Heyer et al. (2009) argued sub-thermal excitation of ^"^CO 
rotational levels may be common in GMCs, causing lower values of line emission than expected under LTE conditions and 
thus causing estimates of E based on LTE assumptions to be underestimates. The IRDC filaments we are considering are 
generally of higher density than the typical GMC volumes considered by Heyer et al. (2009). For the average and highest 
intensity spectra in each IRDC (see Fig. [5]), we evaluate /3 (see Table [T]). We also estimate 71h2, assuming a filament 
line-of-sight thickness that is equal to its observed width. A comparison of nH2 with /?ncrit shows that the densities are 
close to or greater than the effective critical densities, thus justifying the assumption of LTE conditions. Furthermore, we 
expect that the typical density at a given location is greater than our estimated values due to clumping on angular scales 
that are smaller than the 46" resolution of the CO observations. Such clumping is apparent in the MIR extinction maps 

(Figs. [31 & a. 

3. COMPARISON Si3CO AND SsMp: POSSIBLE TRENDS IN TEMPERATURE, CO DEPLETION AND DUST OPACITY WITH T, 

The morphologies of IRDCs F and H are shown in Figures [3] and [H respectively. The ^"^CO emission is more extended 
than the structure traced by the BT09 MIR extinction maps, which are derived by comparing the observed 8/.tm intensity 
at the cloud position with the expected background intensity interpolated from nearby regions. Uncertainties in this 
estimation of the background lead to a lack of sensitivity of the MIR extinction maps to mass surface density contrasts 
AE < 0.013 gcm-2. 

Since there is a significant amount of molecular gas in the surrounding "envelope" regions around the filaments, to 
compare the ratio of Esmf and S13C0 in the filaments we must first subtract the contribution to S13C0 from the envelopes. 
To do this we consider orthogonal strips across each IRDC filament and measure ^^CO emission on either side (see Figs. [3] 
&[4]). The spectra of these "off-source", envelope regions are corrected for optical depth effects, as described above, 
averaged and then subtracted from the central "on-source", filament region (see Fig. [5]). We find the envelope-subtracted 



Table 1 

Escape Probabilities and Critical Densities 



IRDC Spectrum A^isco S13CO /3 nH2 /S^crit "H2/{/3ncrit) 

{seeFig.O {10^^ cm-^) (gcm-^) (cm-^) (cm-^) 



Fmean 2.73 0.0344 0.873 1340 1660 0.807 

Fhigh 4.94 0.0622 0.705 2430 1340 1.81 

Hmcan 2.75 0.0347 0.696 2850 1322 2.15 

Hhigh 3.46 0.0436 0.587 3590 1120 3.22 
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Fig. 2. — Column density profiles with velocity for IRDCs F (left) and H (right). The dashed lines shows the profiles assuming optically 
thin emission, while the solid lines show the profiles after correction for optical depth. A temperature of 20 K was assumed for IRDC F and 
15 K for IRDC H (see text). For each cloud the upper, offset profile is that of the average for the cloud and the lower profile is that of the 
highest column density position. 



spectra are generally narrower than the total. We also find that the size of negative residuals created in the subtraction 
process are relatively small. However, the fact that in some cases we see quite significant variation in the envelope spectra 
from one side of the filament to the other, suggests that this is one of the major sources of uncertainty in measuring the 
molecular emission properties of these embedded filaments. After this procedure, we are now in a position to compare 
the envelope-subtracted values of Sisco with those derived from the SMF extinction mapping method of BT09. 

We also note that the MIR extinction derived values of E suffer from their own systematic uncertainties, including 
corrections due to foreground dust emission (BT09), scattering in the IRAC array (Battersby et al. 2010), adopted dust 
opacities and dust to gas ratios. However, saturation due to large 8/im optical depth and/or an insufficiently subtracted 
foreground is not an important source of error for our present study, since the values of Ssmf are in any case reduced 
when the maps are smoothed to the lower resolution of the ^^CO data. 

As noted by BT09, the MIR extinction mapping technique fails for locations where there are bright MIR sources. If the 
intensity of the source is greater than the background model, then formally an unphysical, negative value of S is returned 
by this method. In the analysis of BT09, negative values of S are allowed up to a certain threshold value to account for 
noise-like, approximately Gaussian, fluctuations in the background intensity. For the more extreme fluctuations caused 
by bright sources, BT09 set E = 0. Thus the effect of a bright MIR region within an IRDC is to cause the extinction 
mapping method to underestimate the true mass surface density. For most point-like MIR sources, this effect is quite 
minor after the extinction maps are averaged to the 22.14" pixel scale of the CO observations. However, in IRDC F we 
identify two regions (indicated in Fig. which are significantly affected by bright MIR sources and exclude them from 
the subsequent analysis. 

The values of E13C0 (after envelope subtraction) and Esmf for each pixel in the filament regions of IRDCs F and H are 
compared in Figure [Bl We note that both these measures of E are independent of the distance to the cloud. Considering 
just the data with E13C0 and Esmf > 0.01 gcm~^, the best fit power law relation Ei3co/gcm~^ = A(EsMF/g cm~^)" 
has a = 0.77 ± 0.18, 0.32 ± 0.28 and A = 0.58 ± 0.44, 0.074 ± 0.069 for IRDCs F and H, respectively. 

Over the range 0.01 < EsMF/gcm^^ < 0.07 wc find that E13C0 — Esmf to within a factor of ~ 2 for the mean values 
in each IRDC filament. The dispersion within an individual IRDC from pixel to pixel is also at about this level. The 
systematic offsets may reflect real systematic variations of the assumed temperature, ^'^CO abundance, dust opacities or 
envelope subtraction method for each IRDC compared to our adopted values and methods. The dispersion may reflect 
local systematic variations and errors introduced by measurement noise. 

In both IRDCs the ratio Ei3co/Esmf decreases with increasing mass surface density as measured by Esmf- We gauge 
the significance of this trend by noting that the above values of a for IRDCs F and H differ from unity by 1.3 and 2.4 
standard deviations, respectively, assuming errors are distributed normally. 

There are several physical processes that may be causing such a trend. A systematic temperature decrease from the 
lower column density, outer regions of the filaments to the higher column density centers would lead us to systematically 
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Fig. 3. — Morphology of IRDC F. Top left: Spitzer GLIMPSE IRAC 8 /xm image, with linear intensity scale in MJy sr~^. The horizontal 
line shows a scale of 3'. The image has 1.2" pixels and the PSF has a FWHM of 2"' Top middle: Mass surface density, SsmFi with linear 
intensity scale in g cm~^, derived from the previous image using the small median filter MIR extinction mapping method of Butler & Tan 
(2009). The inner, red rectangle (centered at I = 34.483°, b = 0.276°, with P.A.= +60° and size 0.0582° by 0.198°) along the filament shows 
the "on source" region we consider to contain the main filamentary structure of the IRDC. The outer, blue rectangle extends to "off source" 
regions we consider to be representative of the surrounding GMC envelope. These rectangles are divided into 7 orthogonal strips to aid in 
the separation of components of CO emission from the filament and GMC envelope. Top right: The same extinction map convolved with a 
Gaussian of 46" FWHM to match the resolution of the CO maps and pixelated to 22"on the same grid as the GRS survey image. The two 
black hatched squares show regions with unreliable measures of SsMF because of the presence of bright MIR sources. These are excluded from 
the comparison with Si3co- Bottom left: Integrated intensity map of ^^CO(l-O) emission over the full velocity range of —5 — 135 km s~^ 
of the GRS survey, with linear intensity scale in Kkms~^. Bottom middle: Integrated intensity map of ^-^00 over the velocity range of 
48 — 65 km s~^, i.e. the gas we believe is associated with the IRDC, with linear intensity scale in K km s~^. Bottom right: Mass surface 
density of the filament derived from ^^CO emission, SiacOi with linear intensity scale in gcm~^. 



underestimate Siaco in the latter, while having little direct effect on dust opacities (although there may be an indirect 
effect via formation of ice mantles, see below). For the temperature ranges we expect to be present, i.e. from ^10-20 K 
(e.g. Pillai et al. 2006 - although note these are based on measurements of NH3, which may trace different conditions to 
those of the ^^CO), this effect becomes important for values of Tb,^ > 4 K (see Fig. [TJj). We illustrate the size of this 
effect for the highest EgMP positions in IRDCs F and H in Fig. [51 A systematic temperature decrease of 5 K in the high 
E regions could remove much of the observed trend. 

Another possibility is that our corrections for the optical depth of the ^'^CO emission are systematically underestimated 
at the higher column density positions. This would be expected if a significant amount of mass is contained in unresolved 
dense cores. The largest optical depth corrections in the highest column density positions presently lead to an increase 
in the estimated column by factors < 2. Future higher angular resolution ^^CO and C^^O observations of these filaments 
are required to investigate this issue further. 

Depletion of CO molecules onto dust grains is known to occur in cold, high volume density gas (Caselli et al. 1999). 
This process could systematically reduce the gas phase CO abundance in the high S regions of the IRDCs. For depletion 
to be fully responsible for the observed trends in IRDCs F and H would require about a factor of 2 depletion as Ssmf 
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Fig. 4.— Morphology of IRDC H. Top left: Spitzer GLIMPSE IRAC 8 /im image, with hnear intensity scale in MJy sr~^ The 
horizontal line shows a scale of 3'. The image has 1.2" pixels and the PSF has a FWHM of 2"' Top middle: Mass surface density, 
SsMF, with linear intensity scale in gcm~'^, derived from the previous image using the small median filter MIR extinction mapping 
method of Butler & Tan (2009). The inner, red rectangle (centered at I = 35.517°, b = -0.275°, with RA.= +62.84° and size 
0.0307° by 0.0637°) along the filament shows the "on source" region we consider to contain the main filamentary structure of the 
IRDC. The outer, blue rectangle extends to "off source" regions we consider to be representative of the surrounding GMC envelope. 
These rectangles are divided into 4 orthogonal strips to aid in the separation of components of CO emission from the filament and 
GMC envelope. Top right: The same extinction map convolved with a Gaussian of 46" FWHM to match the resolution of the CO 
maps and pixelated to 22"on the same grid as the GRS survey image. Bottom left: Integrated intensity map of ^^CO over the full 
velocity rang e of -5 - 135 km s"^ of the GRS survey, with linear intensity scale in K km s ^ . Bottom middle: Integrated intensity 
map of "CO over the velocity range of 40 — 50 km s ^, i.e. the gas we believe is associated with the IRDC, with linear intensity 
scale in Kkms~^. Bottom right: Mass surface density of the filament derived from ^"'CO emission, Eiaco, with linear intensity 
scale in g cm~^. 

increases from 0.01 to 0.05 gcm^^. Because of depletion, CO is not expected to be an ideal tracer of the densest, coldest 
parts of IRDCs. High resolution studies of other tracers, such as NH3, will likely be needed to measure the kinematics of 
these regions. 

The depletion of CO via formation of CO ice mantles on dust grains would also have some effect on the MIR opacities 
of these grains, thus affecting our measurement of Ssmf. The grains would become larger and absorption features due 
to pre-existing water ice mantles may become obscured. The Ossenkopf & Henning (1994) grain models show a ~50% 
increase in 8/Lim opacity, KsAim, going from bare grains to those with thick ice mantles. The BT09 estimates of Ssmf 
assumed a constant value of Ks/^m consistent with the thin ice mantle model of Ossenkopf & Henning (1994). Thus the 
observed decrease in the ratio of Eisco/^smf with increasing mass surface density could be caused by thickening of grain 
ice mantles, causing us to systematically overestimate Ssmf- 

Another possible explanation for a trend of decreasing Sisco/^smf with increasing mass surface density is active 
chemical fractionation (Langer et al. 1980, 1984; Glassgold, Huggins, & Langer 1985; Visser, Van Dishoeck, & Black 
2009), which enhances the abundance of ^^CO in regions where ^'^0+ is present via the via the ion-molecule exchange 




Fig. 5. — Subtraction of IRDC envelopes, (a) Top left: Filament and envelope of IRDC F, with the 7 sets of spectra (top to bottom) 
corresponding to the 7 orthogonal strips (top left to bottom right) shown in Figure O In each, the dotted, red line shows the total ^"^CO 
column density distribution, including optical depth corrections, from the filament region, which includes an assumed contribution from 
"envelope" material along the line of sight. The long dashed and dot-dashed, blue lines show the total ^''CO column density from the 
northern and southern envelope regions, respectively, (b) Top right: Filament and envelope of IRDC H, with the 4 sets of spectra (top to 
bottom) corresponding to the 4 orthogonal strips (top left to bottom right) shown in Figure |4] The line styles have the same meaning as 
in (a), (c) Bottom left: For the same strips as in (a), we subtract the average of the northern and southern envelope spectra (short dashed 
blue lines) from the filament (dotted red line), to leave an estimate of the material in the filament (solid, black line). The vertical solid line 
indicates the mean velocity, (d) Bottom right: Same analysis and labels as (c) applied to IRDC H. 
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Fig. 6. — (a) Top Left: Direct comparison of E13C0 (after envelope subtraction) and Esmf in IRDC F. The crosses show locations 
where both Eisco and Esmf > 0.01 g cm~^. The dots show locations of lower surface densities. The uncertainties in the individual 
measurements are assumed to be 15% plus a systematic error of 0.01 gcm~^. The dotted line shows the one to one linear relation 
and the long dashed line shows the best-fit offset linear relation. The solid line shows the best-fit power law relation (see text). 
(b) Top Right: Same as (a) but for IRDC H. (c) Bottom Left: Logarithm of the ratio of E13C0 to Esmf as a function of Esmf for 
IRDC F, with the same symbol notation as in (a). The cross in the upper-right corner indicates typical estimated uncertainties. 
The solid (dashed) arrow shows the effect on the highest E position of reducing the assumed temperature from 20 K to 15 K 
(10 K). (d) Bottom Right: Same as (c) but for IRDC H. The arrow shows the effect on the highest E position of reducing the 
assumed temperature from 15 K to 10 K. 
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reaction (Watson, Anicich, & Huntress 1976), ^^CO + ^^C+ ^-^0+ + ^^CO + 35K. FUV irradiation maintains a relatively 
high abundance of ^^0+ in the outer regions of the cloud, with Ay < 1 mag. The relative abundance of ^^CO to ^^CO 
can be enhanced by factors of ^ 10. For Ay > 2, the enhancement factor is only about 20% (e.g. for a model with 
riHa — 10'^ cm^'^ and T=20 K; Glassgold et al. 1985). For IRDCs F and H, the filaments are typically embedded inside 
"envelope" gas with Sisco — 0.01 g cm~^, i.e. A^h — 4 x 10^^ cm^^. The trend of decreasing Si3co/5^Smf is observed to 
occur as Ssmf increases to ^ 0.05g cm~^. ^^CO fractionation enhancements should not be varying in this regime by large 
enough factors to explain our observed results. It should also be noted that even in the lower-column density regions, if 
the FUV flux is large enough, then isotope selective photodestruction of ^"^CO compared to the self-shielded ^^CO (Bally 
& Langer 1982) may reverse the fractionation-produced enhancement of the abundance of ^'^CO to ^^CO. 

4. THE DYNAMICAL STATE OF THE IRDC FILAMENTS 

4.1. IRDC Masses from "'^^CO Emission and MIR Dust Absorption 

We calculate IRDC masses using the observed mass surface densities and angular sizes, and by assuming near kinematic 
distances, since the clouds are seen in absorption against the Galaxy's diffuse MIR emission. We adopt the kinematic 
distances of Simon et al. (2006), who assumed the Clemens (1985) rotation curve. This leads to a distance of 3.7 kpc for 
IRDC F and 2.9 kpc for IRDC H. We assume uncertainties of 0.5 kpc, which could result, for example, from line of sight 
noncircular motions of '-^ 8 kms~^. Temperature uncertainties of 5 K would lead to Siaco uncertainties of 20%. We 
estimate similar levels of uncertainty in Ssmf due to foreground correction and background interpolation uncertainties 
(BT09). Then, for IRDC F, the ^'^CO-derived mass assuming T = 20 K in the on-source filament region after envelope 
subtraction is Afisco = 3300 ± 1100 M© and the dust extinction mass is A/smf = 1900 ± 640 Mq. For IRDC H, we find 
Mi3co{T = 15 K) = 370 ± 150 M© and Msmf = 580 ± 230 Mq. For our calculations involving the total mass of the 
clouds we take averages of the above estimates while still assuming 20% uncertainties in the averaged values of S. Thus 
we adopt M = 2600 ± 870 M© for IRDC F and M = 480 ± 190 M© for IRDC H. These resuhs are summarized in Tabled 

4.2. Ellipsoidal Cloud Virial Analysis 

Following Bertoldi & McKee (1992, hereafter BM92), we consider an ellipsoidal cloud with radius R normal to the axis 
of symmetry and size 2Z along the axis. The aspect ratio is defined as y = Z/R, while -Rmax and i?min are the semimajor 
and semiminor axes of the ellipse obtained by projecting the cloud onto the plane of the sky. 

IRDCs F and H both have relatively thin, filamentary morphologies: we set i?max/-Rmin — 3.40, 1.78, respectively, i.e. 
the same elongation as the rectangular regions we consider to define the filaments (Figs. Given these morphologies, 

we expect the symmetry axes of the clouds to be close to the plane of the sky. Thus for both we adopt a fiducial value 
of the inclination angle between the cloud symmetry axis and the line of sight of 6* = 60°. We assume R = Rmin and 
Z = i?inax/sin6', so y = 3.93, 2.06 for IRDCs F and H, respectively. An uncertainty of 15° in inchnation would cause 
^ 15% uncertainties in y. It is also useful to introduce a geometric mean observed radius, i?obs = (-Rmax-Rmin)"'^^^, which 
is also related to R via, = Rcos^/^e[l + (ytan6')2]i/4^ bm92 also define Rm as the mean value of i?obs averaged over 
all viewing angles, but for our individual clouds we will express quantities in terms of R and i?obs- It should be noted that 
the treatment of these IRDC filaments as simple ellipsoids is necessarily approximate, and we consider the filamentary 
analysis of section l4?3l to be more accurate. 

Table 2 

IRDC Ellipsoidal Virial Analysis 



Cloud property IRDC F IRDC H 



d (kpc) 


3.7 ±0.5 


2.9 ±0.5 


R (pc) 


1.88 ±0.25 


0.91 ±0.16 


Robs (pc) 


3.52 ±0.48 


1.21 ±0.21 


y = Z/R 


3.93 ± 0.6 


2.06 ±0.3 


Si3CO (gem 2) 


0.0147 ±0.003 


0.0133 ±0.003 


SSMF (gcm~2) 


0.0082 ± 0.0016 


0.0209 ± 0.004 


Si3Co(env) (gcm-2) 


0.0196 ± 0.004 


0.0212 ±0.004 


Mi3co (Me) 


3300 ±1100 


370 ± 150 


Msmf (Mq) 


1900 ± 640 


580 ± 230 


M (Mq) 


2600 ± 870 


480 ± 190 


ai 


10/9 


10/9 


02 


0.538 


0.750 


W (lO^i^erg) 


-11.0 ±7.1 


-1.08 ±0.78 


cr (km/s) 


1.46 ±0.15 


1.20 ±0.12 


ts = 2R/cT (Myr) 


2.5 ±0.4 


1.5 ±0.3 


r (lO-^^erg) 


16.4 ±6.4 


2.06 ±0.91 


ToiA) (lO^erg) 


14 ±8 


2.7± 1.8 


To{B) (lO^erg) 


73 ±30 


4.3 ±2.0 


a = bcr'^R/{GM) 


1.78 ±0.73 


3.17 ± 1.50 
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If the clump is in an environment that is evolving with a time scale longer than the clump's dynamical time scale or 
signal crossing time tg = 2R/a, then it should obey the equilibrium virial equation (McKee & Zweibel 1992): 

= 2{T-To) + M + W. (8) 

Here, T is the clump kinetic energy, To is the kinetic energy resulting from the surface pressure on the clump, A4 is the 
magnetic energy associated with the cloud, and W is the gravitational binding energy, which for an ellipsoidal cloud is 
(BM92) 

3 GM2 

W = --aia2^-, (9) 

where for a power-law density distribution p oc r~'''' , ai = (1 — kp/3)/{l — 2kp/5) and 

arcsinh(y2 _ 1)1/2 

«2 = _ 1)1/2 (10) 

for prolate clouds. Note, our definition of W and 02 differs slightly from BM92 since we do not need to consider Rm- 
We adopt kp = 1 (based on a study of the density profiles in IRDCs - Butler & Tan, in prep.) so that ai = 10/9. For 
our measured values of y, we have 02 = 0.53, 0.71 for IRDCs F and H, respectively. For the mass of the cloud we take 
the average values of the estimates from MIR dust extinction and ^^CO line emission. Using these values we estimate 
W = —(11.0, 1.08) X 10**^ erg for IRDCs F and H, respectively (see Table 2). The uncertainties in W are relatively large 
given the measurement errors of S and R. 

The clump kinetic energy is 7" = (3/2)Mcr^, where a is the average total ID velocity dispersion, which we derive from 
the "^^CO line emission (counting only those parts of the envelope-subtracted spectra with positive signal greater than or 
equal to one standard deviation of the noise level) including corrections for the molecular weight of ^'^CO and our adopted 
cloud temperatures. We estimate that we measure cr to a 10% accuracy. We find T = (16.4, 2.06) x 10'*^ erg for IRDCs F 
and H, with uncertainties at about the 40% level. 

The surface term for the kinetic energy is To = {3/2)PoV. We estimate this term in two ways. First (method A), we 
can measure the mass surface density of the surrounding molecular cloud from the ^^CO emission of the envelope regions. 
We find Ei3co(env) = 0.0196, 0.0212 g cm'^ for IRDCs F and H, respectively. We scale these values by M/Muco, 
i.e. 0.79 and 1.30. If the envelope is self-gravitating, it has mean internal pressure P(env) = 1.85GS^(env), adapting 
the analysis of McKee & Tan (2003) with parameters fg = l (i.e. fully gas dominated), 0geom = ^obs/(^^^) ~ ^-^ 
(adopting an intermediate value for the two IRDCs accurate to about 20%), = 2.8 (the fiducial value of McKee & 
Tan, measuring the ratio of the total pressure including magnetic fields to that assuming they were absent) and avir = 1- 
Thus setting Pq = P(env) = (2.96,9.38) x 10"^^ cgs for IRDCs F and H and with the cloud volume V = A-kR^Z/2,, we 
find %{A) = (14,2.7) X lO'"^ erg. 

Second (method B), we estimate the density in the envelope region, assuming it has a cylindrical, annular volume 
with outer radius 2R and inner radius R. For IRDCs F and H, we find densities of p = 4S(env)/(37ri?) = (1.43,3.20) x 
10~^^ gcm^'*, equivalent to nH(env) = (610, 1600) cm~^. We again scale these values by M/Mi^co, i-e. 0.79 and 1.30 
for IRDCs F and H, respectively. We then equate Po = /3(env)(7^(env), where (j(env) is the velocity dispersion of the 
envelope gas (we find 3.65, 1.89 km s^^ for IRDCs F and H) and evaluate %{B) = {?>/2)VPq (73,4.3) x lO'"^ erg, with 
V = (4/3)7rp2^. 

These results indicate that, for both IRDCs, the surface pressure term of the virial equation is comparable to or much 
larger than the internal kinetic term, although the uncertainties are large. Assuming To > T, then for virial equilbrium 
to be maintained would require AA = (l/87r) jy {B^ — Bf^dV > —W, where Bq is the magnetic field strength far from the 
cloud and Va is a volume that extends beyond the cloud where the field lines have been distorted by the cloud. Assuming 
Va is the volume of the envelope regions and assuming negligible Bq, we find B > 10, 13 /uG for IRDCs F and H. If a 
more realistic value of Bq = 10 /iG is adopted (Crutchcr et al. 2010), then we find B > 14, 16 /xG. Thus relatively modest 
magnetic field enhancements could stabilize the clouds. 

Bertoldi & McKee (1992) define a dimensionless virial parameter, a = 5(j'^ Rm/ (GM) to describe the dynamical state 
of clouds. We adopt a slightly revised definition a = 5a'^R/{GM) = aia22T/\W\, with 02 defined as above. For IRDCs F 
and H we find a = 1.78, 3.17. While these values are quite close to unity, especially for IRDC F, which is commonly taken 
to infer that self-gravity is important, this is somewhat misleading since the value of 02 is quite small for these elongated 
clouds and the surface pressure terms seem to be quite important. 

It is possible that these IRDCs have not yet reached virial equilibrium if their surroundings are evolving on timescales 
< is ^ a few Myr. The mean velocity of the ^"^CO emitting gas in the north and south envelope regions is 55.72, 
56.10 km s""'^, respectively, for IRDC F, and 44.96, 44.28 km s""'^, respectively, for IRDC H. The north/south velocity 
dispersions are 3.81/3.38 km s"^ for F and 1.69/2.07 km s"^ for H. The north/south Siaco's arc 0.040/0.031 gem for 
F and 0.050/0.051 g cm~^ for H. The ratios of north/south pressures (estimated by method B, oc Ecr^) are thus 1.64 and 
0.65 for IRDCs F and H. The pressures appear to be fairly similar on the different sides of the filaments, although the 
uncertainties are such that variation at the level of about a factor of two could be present. 
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4.3. Filamentary Cloud Virial Analysis 

Fiege & Pudritz (2000, hereafter FPOO) present a virial analysis of filamentary clouds. They derived the following 
equation satisfied by pressure-confined, nonrotating, self-gravitating, filamentary (i.e. lengths 3> widths) clouds threaded 
by helical magnetic fields that are in virial equilibrium; 

(11) 



Here Pq is the external pressure at the surface of the filament, P = pa^ is the average total pressure in the filament, m 
is the mass per unit length, rrivir = /G is the virial mass per unit length. Mi is the magnetic energy per unit length, 
and Wi — ^m?G is the gravitational energy per unit length. 

We divide IRDCs F, H into 7, 4 orthogonal strips (see Figs.|n]&lll) with angular widths 1.70', 0.955' along the filaments, 
respectively. Assuming a fiducial value of the inclination angle between the cloud symmetry axis and the line of sight of 
6 = 60°, these correspond to physical lengths along the filaments of 2.11, 0.930 pc, respectively. 

In Table[3l for each strip in IRDCs F and H, we list the values of E13CO1 SSMF7 M (calculated from the mean of these 
values of E), to, p, v, a, TOvh-, P, Si3co(env), p(env) (calculated after scaling I]i3co(env) by 0.5(Si3co + Esmf)/Ei3co), 
z;(env), (T(env) and P(env) = p(env)(T^(env). We equate Pq = P(env). 

Following FPOO, in Fig. [7] we plot Pq/P versus m/ rrivir- The range of models considered by FPOO allows for positive 
values of M-i/\Wi\ (i.e. poloidally-dominated B-fields that provide net support to the filament against gravitational 
collapse) and negative values (i.e. toroidally-dominated B-fields that provide net confinement of the filament). In all 
cases, Pq/P < 1. In contrast, we find all of the filament regions have Pq/P > 1, i.e. the pressures in the envelope regions 
appear to be greater than in the filament. This echoes the results from the ellipsoidal virial analysis, which found large 
surface pressure terms. Assuming our measurements of pressures are reliable, e.g. are not being adversely affected by 
systematic effects due to, for example, our assumed filament and envelope geometry, then these results imply that the 
filaments have not yet reached virial equilibrium. 

Very large values of Pq/P are inferred for strips Fl, F2 and F3. These are consistent with the filament and envelope 
spectra shown in Fig.[5l which reveal a relatively weak filament and relatively strong and varying envelope velocity profiles. 

Strips F6, F7 and 112 have the smallest values of Pq/P < 2. Examining the IRAC 8/im images (Figs. O&S]), there 
is some indication that these are the sites of relatively active star formation (especially F7 and 112). Star formation 
requires gravitationally unstable conditions in the filament, i.e. regions where self-gravity starts to dominate over external 
pressure. Our results indicate that this also requires the local region of the filament to reach approximate virial equilibrium, 
although surface pressure terms still remain dynamically important, i.e. m/mvir is significantly less than unity. Given our 
measurement uncertainties and the fact that the observed regions do not have Pq/P < 1, we are not able to determine 
whether the field geometries are more dominated by poloidal or toroidal components. We note that FPOO's conclusion 
that observed filaments are dominated by toroidal fields depends on their assumption that the filaments are in virial 
equilibrium and on their choice of Pq, which was not directly measured for most of the sources they considered. 

We caution that if independent molecular clouds are present along these lines of sight and with similar velocities to the 
filaments, then this may cause us to overestimate the velocity dispersion and pressure in the envelope regions around the 
filaments. From the spectra shown in Fig. [5] we do not expect this is occurring in IRDC H, since the envelope spectra 
share a very similar velocity range as the filament. The situation in IRDC F is less clear cut, since there appears to be a 
broader, offset component that dominates more in the envelope region. 

5. CONCLUSIONS 

We have compared measurements of mass surface density, E, in two IRDC filaments based on ^"^CO observations, E13C0, 
with those derived from MIR extinction mapping, Ssmf, finding agreement at the factor of ^ 2 level. A systematic decrease 



Table 3 

IRDC Filamentary Virial Analysis 



Cloud property 


Fl 


F2 


F3 


F4 


F5 


F6 


F7 


Ftot 


HI 


H2 


H3 


H4 


Htot 


Sl3CO (10-2gcm-2) 


0.218 


0.338 


0.811 


1.31 


2.29 


2.87 


2.58 


1.47 


0.864 


1.67 


1.56 


1.21 


1.33 


SsMF (10-2gcm-2) 


0.476 


0.311 


-0.0432 


0.862 


1.37 


1.62 


1.03 


0.823 


1.38 


1.94 


1.97 


3.28 


2.09 


M (Mq) 


114 


107 


126 


357 


601 


737 


593 


2640 


78.5 


126 


124 


157 


478 


m (M0pc~l) 


54.0 


50.7 


59.7 


169 


285 


349 


281 


178 


84.4 


135 


133 


169 


128 


p (10~^-^g cm~^) 


3.29 


3.09 


3.64 


10.3 


17.4 


21.3 


17.1 


10.9 


22.0 


35.1 


34.6 


44.0 


33.3 


V (km s"-*^) 


59.19 


59.34 


58.81 


58.26 


58.30 


57.77 


57.77 


58.39 


45.09 


45.20 


45.13 


44.78 


45.07 


a (km s~^) 


1.36 


1.27 


1.07 


1.44 


1.59 


1.86 


2.27 


1.46 


1.52 


1.34 


1.03 


0.995 


1.20 


mvir (M0pc~i) 


856 


751 


536 


962 


1170 


1610 


2400 


986 


1070 


840 


494 


460 


669 


P (10-i2cgs) 


6.05 


4.99 


4.19 


21.3 


43.8 


73.8 


88.3 


23.1 


50.7 


63.4 


36.8 


43.6 


47.9 


Sl3Co(env) (10-2gcm-2) 


1.61 


1.83 


1.94 


2.14 


1.81 


1.95 


2.46 


1.96 


2.31 


1.86 


2.00 


2.35 


2.12 


p(env) (10~^^ g cm~'^) 


18.7 


12.9 


6.72 


13.0 


10.6 


11.2 


12.6 


11.2 


45.3 


30.4 


34.2 


65.9 


41.2 


v{env) (km s"-') 


56.03 


55.84 


55.78 


56.59 


56.04 


55.65 


55.72 


55.95 


44.97 


44.78 


44.53 


44.40 


44.67 


cr(env) (km s"-') 


3.45 


3.48 


3.95 


3.91 


3.55 


3.43 


3.56 


3.65 


1.85 


1.91 


1.89 


1.86 


1.89 


P(env) (lO-l^cgs) 


223 


156 


105 


199 


133 


132 


159 


149 


155 


110 


122 


229 


148 
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Fig. 7. — Po/P versus m/m^i^ for strips in IRDCs F (open squares joined by dotted li ne f or Fl to F7) and H (open triangles joined by 
dashed line for HI to H4). The smooth curves show the conditions satisfied by equation for Mi/\Wi\ < (solid lines), Mi/\Wi\ = 
(dotted line), and A^;/|W;| > (dashed lines). The large observed values of Po/P may indicate that most of the regions in IRDCs F and H 
have not yet established virial equilibrium. Those strips v^ith the lowest values of Pq/P, i.e. F6, F7, and H2, appear to be undergoing more 
active star formation (see text and Figs. [3l&:l4)l. 
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in Siaco/^SMF with increasing E may be due to a systematic decrease in temperature, increase in the contribution of 
under-resolved high optical depth regions, increase in dust opacity, or decrease in ^■^CO abundance due to depletion in 
regions of higher column density. Future studies that spatially resolve the temperature structure and MIR dust absorption 
properties can help to distinguish these possibilities. 

We have then used the kinematic information derived from ^'^CO to study the dynamical state of the IRDCs. In 
particular we have evaluated the terms of the steady-state virial equation, including surface terms, under the assumption 
of ellipsoidal and filamentary geometries. In both cases we find evidence that the surface pressure terms are important 
and possibly dominant, which may indicate that the filaments, at least globally, have not yet reached virial equilibrium. 

These results would be consistent with models of compression of dense gas in colliding molecular flows, e.g. GMC 
collisions. Tan (2000) proposed that this mechanism may trigger the majority of star formation in shearing disk galaxies. 
The expected collision velocities are ^ lOkms"^. It is less clear whether colliding atomic flows, (e.g. Heitsch et al. 2008), 
which form the molecular gas after shock compression of atomic gas, would also produce such kinematic signatures: recall 
that we are inferring large surface pressures based on ^'^CO emission from the envelopes around the IRDC filaments. 

Recent observations of extended, parsec-scale SiO emission, likely produced in shocks with velocities > 12 kms~^ in 
IRDC H by Jimenez-Serra et al. (2010) may also support models of filament formation from converging fiows. However, 
we caution that the observed extended SiO emission is very weak and may also be produced by multiple protostellar 
outflow sources forming in the IRDC (see Jimenez-Serra et al. 2010 for further discussion). 

Our resolved fllamentary virial analysis also indicates that the regions closest to virial equilibrium (strips F6, F7 and 
H2) are those which have initiated the most active star formation. This would be expected if models of slow, equilibrium 
star formation (Tan et al. 2006; Krumholz & Tan 2007) apply locally in these regions. In this case, these dense regions 
that have become gravitationally unstable, perhaps due to the action of external pressure and/or converging flows, then 
persist for more than one local dynamical time and so are able to reach approximate pressure and virial equilibrium with 
their surroundings. In this scenario, they are stabilized by the ram pressure generated by protostellar outflow feedback 
from the forming stars (Nakamura & Li 2007). 

We thank Michael Butler for providing the MIR extinction maps used in this analysis. We also thank Peter Barnes, 
Crystal Brogan, Michael Butler, Paola Caselli, James Jackson, Izaskun Jimenez-Serra, Thushara Filial and Robert Simon 
for helpful discussions. The comments of an anonymous referee led to improvements in the paper. AKH acknowleges 
support from a SEAGEP Dissertation Fellowship. JCT acknowledges support from NSF CAREER grant AST-0645412 
and NASA Astrophysics Theory and Fundamental Physics grant ATP09-0094. 

REFERENCES 



Battersby, C, Bally, J., Jackson, J. M. et al. 2010, ApJ, 721, 222 

Bally J., & Langer, W. D. 1982, ApJ, 255, 153 

Bertoldi, P., & McKee, C. F. 1992, ApJ, 395, 140 {BM92) 

Butler M. J., & Tan J. C. 2009, ApJ, 696, 484 (BT09) 

Carey S. J., Clark, F. O., Egan, M. P., Price, S. D., Shipman, R. F., 

& Kuchar, T. A. 1998, ApJ, 508, 721 
Carey S. J., Feldman, P. A., Redman, R. O., Egan, M. P., MacLeod, 

J. M., & Price, S. D. 2000, ApJ, 508, 721 
Caselli, P., Walmsley C. M., Tafalla, M., Dore, L., & Myers, P. C. 

1999, ApJ, 523, L165 
Clemens, D. P. 1985, ApJ, 295, 422 

Crutcher, R. M., Wandelt, B, & Falgarone, E. 2010, in From 
Stars to Galaxies, J. Tan & S. Van Loo (eds.), e-proceedings: 
|http://conference^astro.ufl.edu/STARSTOGALAXIES/ p3 
Devme, K. (ZU09), f h.U. thesis, University of Wisconsin 
Fiege, J. D., & Pudritz, R. E. 2000, MNRAS, 311, 85 (FPOO) 
Glassgold, A. E., Huggins, P. J., & Langer, W. D. 1985, ApJ, 290, 
615 

Goodman A. A., Pineda J. E., Schnee, S. L. 2009, ApJ, 692, 91 
Heitsch, F., Hartmann, L. W., Slyz, A. D., Devriendt, J. E. G., & 

Burkert, A. 2008, ApJ, 674, 316 
Heyer, M. H., Krawczyk, C, Duval, J., & Jackson, J. M. 2009, ApJ, 

699, 1092 

Jackson, J. M., Rathborne, J. M., Shah, R. Y. et al. 2006, ApJS, 
163, 145 

Jimenez-Serra, I., Caselli, P., Tan, J. C, Hernandez, A. K., Fontani, 
F., Butler, M. J., & Van Loo, S. 2010, MNRAS, 406, 187 



Krumholz, M. R., & Tan, J. C. 2007, ApJ, 654, 304 
Lacy, J. H.; Knacke, R.; Geballe, T. R.; Tokunaga, A. T. 1994, ApJ, 
428, L69 

Langer, W. D., Goldsmith, P. F., Carlson, E. R., & Wilson, R.W. 

1980, ApJ, 235, L39 
Langer, W. D., Graedel, T. E., Frerking, M. A., & Armentrout, P. 

B. 1984, ApJ, 277, 581 
McKee, C. F., & Tan, J. C. 2003, ApJ, 585, 850 
McKee, C. F., & Zweibel, E. G. 1992, ApJ, 399, 551 
Milam, S. N., Savage, C, Brewster, M. A., & Ziurys, L. M., 2005, 

ApJ, 634, 1126 
Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 
Nakamura, F. & Li, Z.-Y. 2007, ApJ, 662, 395 

Pillai, T., Wyrowski, F., Carey S. J., & Menten 2006, A&A, 450, 
569 

Pineda, J. E., CaseUi, P., & Goodman, A. A. 2008, ApJ, 679, 481 
Rathborne, J. M., Jackson, J. M., & Simon, R. 2006, ApJ, 641, 389 
Simon, R., Rathborne, J. M., Shah, R. Y., Jackson, J. M., & 

Chambers, E. T. 2006, ApJ, 653, 1325 
Tan, J. C. 2000, ApJ, 536, 173 

Tan, J. C. 2007, in Triggered Star Formation in a Turbulent ISM, 

eds. B. G. Elmegreen & J. Palous, lAUS 237, CUP, p258 
Tan, J. C, Krumholz, M. R., McKee, C. F. 2006, ApJ, 641, L121 
Visser, R., van Dishoeck, E. F., & Black, J. H., 2009, A&A, 503, 323 
Watson, W. D., Anicich, V. G., & Huntress, W. T., Jr. 1976, 205, 
L165 



